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Abstract 

Kerr oscillators are model systems which have practical applica- 
tions in nonlinear optics. Optical Kerr effect i.e. interaction of optical 
waves with nonlinear medium with polarizability x^^^ is the basic phe- 
nomenon needed to explain for example the process of light transmis- 
sion in fibers and optical couplers. In this paper we analyze the two 
Kerr oscillators coupler and we show that there is a possibility to con- 
trol the dynamics of this system, especially by switching its dynamics 
from periodic to chaotic motion and vice versa. Moreover the switch- 
ing between two different stable periodic states is investigated. The 
stability of the system is described by the so-called maps of Lyapunov 
exponents in parametric spaces. Comparison of basins of attractions 
between two Kerr couplers and a single Kerr system is also presented. 



Keywords: Kerr effect, Kerr couplers, Lyapunov exponents, basins of 
attractions, control of system dynamics, chaotic beats. 
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1 Introduction 



One of the best known and most intensively studied optical models is an 
oscillator with Kerr nonlinearity. Different kinds of anharmonic Kerr oscil- 
lators have also been used to study classical and quantum chaos fT]- [5]. 
Mutually coupled Kerr oscillators can be successfully used for a study of 
couplers, the systems consisting of a pair of coupled Kerr fibres. The first 
two-mode Kerr coupler has been proposed by Jensen [6] and investigated in 
depth in [SI E] • Kerr couplers affected by quantization can exhibit various 
quantum properties such as squeezing of vacuum fluctuations, sub-Poissonian 
statistics, collapses and revivals [8l|9]. 

In the last two decades since the publication of the paper by Pecora and 
Carroll [10] the phenomenon of synchronization in systems of the coupled os- 
cillators has become a subject of comprehensive investigation. The problem 
of synchronization of two linearly coupled Kerr oscillators has been studied 
in [11] and the possibility of synchronization of chaotic motion was proved 
numerically. Moreover the case of synchronization of two kinds of Kerr cou- 
plers having a structure of low-dimensional chains (ring and open) has been 
analysed [T2| . 

This paper is an attempt at using the modern tools of nonlinear science 
for numerical investigation of dynamics of a system made of two coupled 
Kerr oscillators. 

In Sec. 2 the basic equation of motion for the single Kerr oscillator is 
introduced. Simple periodic solutions of equations of motions have been 
found and the dynamics of the system as well as basins of attraction for 
such solutions are investigated. Moreover, we calculate the Lyapunov maps 
for the single Kerr oscillator with external periodic as well as modulated 
fields. The second case of the external field is used to generate the so-called 
chaotic beats. In Sec. 3 our single Kerr system analysis is extended over 
the case of two coupled Kerr sub-systems with nonlinear coupling. We find 
the analytic periodic solutions of such a system. Lyapunov maps and the 
basins of attraction for this system are helpful tools in analysis of properties 
of these system. We find that it is possible to change the periodic states of 
one sub-system by changing the initial conditions of the other sub-system 
(switching the periodic dynamics). Moreover, it is proved that coupled Kerr 
sub-systems are able to generate of chaotic beats. 

2 The single Kerr oscillator 
2.1 Equations of motion 

We study the dynamical system described by the following hamiltonian: 

H = Ho + H^, (1) 



2 



where: 

Ho = uja*a + ^ea*^a^ , (2) 

Hi = iF (a*e-'^^' - ae'^^') . (3) 

The hamiltonian Hq represents the so called Kerr oscillator (if e = 0, then Hq 
refers to the harmonic oscillator), whereas the hamiltonian Hi describes the 
interaction of the Kerr oscillator with the periodic external field. The quan- 
tities a and a* are complex dynamical variables describing the amplitudes, 
u denotes the frequency of the free vibrations of the harmonic oscillator - 
basic frequency, e is the parameter describing the Kerr nonlinearity in the 
system (this is the nonlinearity of the third order) and F is the external field 
amplitude at the frequency Vtp. 

The equation of motion for variable a has the form: 

^ = -iioa - iea*a^ + Fe'^^"^ - 7a . (4) 

The term —7a - added on phenomenological grounds - describes the mecha- 
nism of loss with the damping constant 7. All the parameters, that is w, e, 
F, Vtp and 7 are taken to be real. The equation of motion for a* is simply a 
complex conjugation of Eq.(jl]). 

In the autonomous and conservative case, that is when 7 = F = 0, the 
solution of equation (jl]) has a well-known form: 

a{t) = aoe-^('^+^''S"o)t ^ (5) 

where a(t) \t=o= Qo is the initial condition. 

In the nonautonomous case of Eq.(j4]) we can find the periodic solution: 

a{t) = a;e-*('"+^"*")* , (6) 

that is in the form of the solution of the autonomous one. Function (|6]) 
satisfies the equation of motion (jl]) provided that x = F/7. As a result the 
periodic solution of equation (j4]) has the form: 



a{t) = -e V . (7) 

7 

Formally, the function in the form of ([7]) is the solution of the differential 
equation (jl]) only if two condition are fulfilled: (A) Qp = u + e^, and 
(B) the initial condition has the form a(0) = ao = F/7. In other words, the 
periodic solution ([7]) is correct only for special choice of the set of parameters 
u, Qp, F, 7. Finally, it is worth noting that in the nonautonomous case the 
period of solution depends on the initial condition, as well as in the 
autonomous case. In the phase plane {Re a, Im a) the periodic solution 
satisfies the phase equation (circle): 

(Rea)^ + {Ima)^ = F^/-f^ . (8) 
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for any values of frequency oj. 

It should be emphasized that the method presented here is usefull to 
find only the one periodic solution for a given set of the system parameters. 
Generally, the Eq. (4) have up to three (not only periodic) solutions. 

2.2 The dynamics of the system in the phase space 

As a numerical example let us consider the dynamics of a system described 
by equation (j4]), if w = 1, 7 = 0.5, e = 0.01, F = 5 and VLp = 2. Then, in 
compliance with ([7]), the periodic solution of equation (jlj) has the form: 

a{t) = lOe-2** , (9) 

and in the phase space it satisfies the following equation: 

{Re af + {Imaf = 100. (10) 

As a result, for the initial condition = 10 the phase point draws simply a 
circle described by equation (ITUl) . But if the system (jl]) starts from another 
initial conditions, we observe the following interesting behaviour: after some 
time the phase point tends to one of the two orbits: {Re a)^ + {Im a)^ = 100 
or {ReaY + {Imaf = 50. Two examples of such behaviour of the phase 
point are illustrated in Fig. [H In Fig. [l](a) the phase point (marked in pink) 
starts from the initial condition Rea = 2, Im a = and after the time t = 50 
goes into the orbit of the radius r = a/50. However in Fig. W(Jo) we can see 
how the phase point starting from the initial condition Rea = 10, Im a = 15 
goes into the orbit of the radius r = 10, described by equation ( fTOl) . It should 
be emphasized that the orbit r = \/50 is also periodic but it does not fulfill 
conditions (A) and (B). 

Both orbits: {ReaY + {ImaY = 100 and {Rea^ + {ImaY = 50 are 
attractors of the system (jlj). It means that for the parameters: u = 1, 
7 = 0.5, e = 0.01, F = 5 and flp = 2 the system (jl]) tends to one of the two 
steady states (periodic), represented by these two orbits. 

It is easy to show that the orbit described by equation (i?e a)^ + {Im a)^ = 
50 is identical with that generated by equation (jl]) if a; = 1, 7 = 0.5, e = 0.01, 
F = 5 and Qp = 1 and for the initial condition: Rea = 5, Im a = 5. This 
solution has the form: a{t) = (5 + 5z)e~^**. 

2.3 Basin of attraction 

To illustrate the full influence of the initial conditions on the evolution of 
the system we used the so called basins of attraction. Basin of attraction 
is the set of initial conditions which lead to the system's attractor. The 
basins of attraction of the system (jl]) for a; = 1, 7 = 0.5, e = 0.01, F = 5 
and Qp = 2 are presented in FigJ21 There are two attractors of the system 
(limit cycles r = 10 and rf = v^50) and their basins of attraction are marked 
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Figure 1: The phase trajectories of the system (jlj) forcj = 1, 7 = 0.5, 
e = 0.01, F = 5 and Qp = 2 and for the initial conditions: (a) Rea = 2, 
Ima = 0; (b) Rea = 10, Ima = 15. 



5 



Im a 



-30.0 



-15.0 



30.0 



15.0 



0.0 




-30.0 -15.0 0.0 15.0 30.0 f^g 



Figure 2: Basins of attraction for two attractors - limit cycles for r = 10 and 
rt = ^/50. The parameters of the system ([Tj) are: a; = 1, 7 = 0.5, e = 0.01, 
F = 5 and n„ = 2. 



by different colours, the yellow area marks the basin of attraction of the 
attractor {Re ai)^ + {Im ai)^ = 100, whereas the blue one refers to the basin 
of attraction of the attractor {ReaiY + {Imaif' = 50. Both basins have 
interesting geometries. The basin corresponding to the circle of the radius 
rf = a/50 has a spiral-like form with the slip and width decreasing when 
moving away from the centre. The remaining area (blue colour) refers to the 
basin of attraction of the second attractor (the circle of the radius r = 10). 
Both attractors have a special property. They are localised in such a way, that 
each attractor is located partly in its own basin of attraction and partly in 
the basin of attraction of the other attractor. As a result - if the phase point 
starts from the part of the attractor situated in the basin of attraction of the 
other one, it escapes to the other attractor. However, if it starts from the 
part of the attractor situated in its own basin of attraction, it does not change 
the attractor. In analogy to semistable orbits these attractors can be called 
the semistable attractors. So, the system with Kerr nonlinearity is tuneable: 
an adequate choice of the initial condition can result in the transition of the 
phase point from the one attractor to the other. This property seems to be 
useful in applications in optical switches. 

2.4 Parameters detuning 

One of the important properties of dynamical systems is their sensitivity to 
the change in the system's parameters. A small change in a parameter can 
lead to radical changes in the dynamics of the system. This feature is fre- 
quently used to control the dynamical systems. Globally, the behaviour of 
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Figure 3: Lyapunov map - the values of the maximal Lyapunow exponent 
Ai for the parameters of the system (11): w = 1, e = 0.01 and for the pump 
field amplitude F = 5. The initial condition is: Rea = 10, Ima = 0. In the 
parameter space {Qp, 7) the colours correspond to appropriate values of Ai. 



the system can be shown on the so-called Lyapunov map in a parametric 
space. We used here the well known procedure [13] for numerical calcula- 
tion of Lyapunov exponents (Lyapunov spectrum Aj,i = 1,2). In Figl3]we 
show the map of maximal Lyapunow exponent Ai for the parameters of the 
system u = 1, e = 0.01 and for the pump field amplitude F = 5. The map 
is presented in the parameters space {flp, 7). The highest values of Ai cor- 
responding to chaotic oscillations are marked by red and blue colours. The 
chaotic motion exists only for weak damping. For higher damping constants 
we find only single islands of chaotic motion in the pump field parametric 
space {flp, F) (see FigH]). 

To control our system of Kerr oscillators we first change the value of one 
of the parameters (7, F or Qp) of the system (jlj) at time ti. In such a way, 
the phase point being in one of the two orbits (attractors) shown in FigJ2] 
escapes from the attractor to the transient state and sometimes to the new 
periodic state. Then after returning at time ^2 > to the initial values of 
this parameter, the additional periodic state disappears and the phase point 
trajectory tends to one of the two attractors of the system, depending on 
which basin of attraction was the phase point at time ^2 ■ So, through the 
appropriate choices of times ti end t2 as well as the values of the parameters 
of the system we can control its evolution. In particular, we can switch the 
system between two stable periodic states. Such situations are illustrated in 
Fig. [S](a)-(c). In FigEl^a) the phase point starting from the orbit of the radius 
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Figure 4: Lyapunov maps - the values of the maximal Lyapunow exponent 
Ai for the parameters of the system (11): u = 1, e = 0.01 and for the initial 
condition: Rea = 10, Im a = 0. The parameter spaces {Qp, F) are presented 
for (a) 7 = 0.1 and (b) 7 = 0.5. For higher values of 7 only the periodic 
behaviour of the system is observed. 



8 



-10 

-10 10 20 -10 -5 5 10 
Re a Re a 




-10 -5 5 10 
Re a 



Figure 5: The phase trajectories of the system ([5]) for the initial conditions: 
Rea = 10, Ima = and for the following values of the system's parameters: 
(a) w = 1, e = 0.01, F = 5, fip = 2 and 7 = 0.5 for < t < 20, 7 = 0.02 
for 20 < t < 40 and 7 = 0.5 for 40 < t < 60; (b) a; = 1, 7 = 0.5, e = 0.01, 
F = 5 and fip = 2 for < t < 20, fip = 4 for 20 < t < 40 and fip = 2 for 
40 < t < 100; (c) w = 1, 7 = 0.5, e = 0.01, = 2 and F = 5 for < t < 20, 
F = 2 for 20 < t < 40 and F = 5 for 40 < t < 250. 
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r = 10 (marked in red) after detuning the value of the damping constant 7 
in time ti = 20 from 7 = 0.5 to 7 = 0.02 escapes through a transient state 
to the new periodic state r = 11, 91. After coming back with 7 to the initial 
value 7 = 0.5 it goes into the initial orbit (of the radius r = 10). In Fig. O^b) 
the phase point starts also from the point lying on the orbit of the radius 
r = 10 and after detuning the parameter Qp from f2p = 2 to f2p = 4 at 
time ti = 20 it escapes through a transient state to the new periodic state 
r = 1, 84. Then after coming back with the parameter Qp at time ^2 = 40 to 
the initial value fip = 2 it goes into the orbit of the radius rf = y/bO - the 
case of switching between the periodic orbits. Moreover, Figj5]^c) shows the 
phase point starting from the site lying on the orbit of the radius r = 10 after 
detuning the value of the parameter F from F = 5 to F = 2 at time ti = 20 
also moves to another periodic orbit r = 1, 84, and then after returning at 
time t2 = 40 to the initial value of parameter F {F = 5) it goes into the orbit 
of the radius rf = a/50. Generally, it is very interesting that if we change the 
parameters of the system, some extra periodic orbits appear. 

2.5 Generation of chaotic beats 

Since the publication of [H], the new type of signals called "chaotic beats" 
has been investigated [I5] and experimentally generated [I6l |T9j. There 
are two basic kinds of chaotic beats: (1) the signals with chaotic envelopes 
and a stable fundamental frequency, and (2) the signals with almost regular 
collapses and revivals with small chaotic perturbations. Generally, when the 
system is subjected to an external field, we can generate the chaotic beats in 
two ways by modulation of the amplitude or frequency of the external field. 

The more effective method of generation of chaotic beats in system (j3]) 
seems to be the frequency modulation, according to the formula: 

Qp — 7- Qp{l + AQp sin(/it)), 

where fi is the frequency of modulation parameter, and AQp is the am- 
plitude of this modulation. Then the equation of motion for variable a takes 
the form: 

^ = -^^a - tea* a' - 7a + Fe-^^''(i+^^-^^'^('^*»* . (11) 

In numerical calculations we put AQp = 0.1. 

The global dynamics of system f|TT]) is presented in Fig. [6] as a Lyapunov 
map in the parametric space of (7, /i), where the values of the first Lyapunov 
exponent Ai are marked by appropriate colours. The highest values of Ai cor- 
responding to chaotic oscillations are marked by red and blue colours. As we 
can see, they are concentrated in the lower part of the map which corresponds 
to the low values of the damping parameter 7 (mainly for 7 < 0.05). For 
higher 7 the map is dominated by the black and grey colours corresponding 
to the periodic states. An example of chaotic beats generated in the system 
( ITT]) through frequency modulation (with 7 = 0.02 and /i = 0.15) is pre- 
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Figure 6: The values of maximal Lyapunov exponent Ai for the system (fTTi) 
as a function of the damping parameter 7 and the frequency modulation /i, 
for e = 0.01, F = 5, Qp = 2 and with the initial condition: Rea = 10, 
Im a = 0. 



sented in FigjTl The spectrum of the Lyapunov exponents {0.0519, —0.0813} 
of that system with the positive value of Ai indicates chaotic behaviour. 

Similar results were obtained for the resonance case {Qp = u). The 
values of maximal Lyapunov exponent Ai for system (fTT!) as a function of 
the damping parameter 7 (0 < 7 < 0.5) and the frequency modulation /x 
(0 < /i < 1) and for Qp = (jj = 1 are presented in FiglHl 

Analogically as for the non-resonance case, the chaotic behaviour is ob- 
tained only for the low values of the damping parameter 7 (mainly for 
7 < 0.05). In Figl9]we can see the time dependence of Rea of the system 
(dH), for 7 = 0.01, /i = 0.15 and Vlp = to = 1. 

The spectrum of Lyapunov exponents of the beats shown in FigjO] is 
{0.1085,-0.1230} and contains positive value of Ai indicating chaotic be- 
haviour. 

3 The two coupled Kerr oscillators 
3.1 Equations of motion 

It is interesting to know what happens to the dynamics of the single Kerr 
oscillator after coupling it nonlinearly to another analogous oscillator but of 
different frequency. Such a system of two nonlinearly coupled Kerr oscillators 
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Figure 7: Time dependence of Re a of the system (fTTl) . for 7 = 0.02 and 
/i = 0.15. The other parameters are: a; = 1, F = 5, e = 0.01 and VLp = 2. 
The system starts from the initial condition: Rea = 10, Ima = 0. Chaotic 
beats. 




Figure 8: The values of the maximal Lyapunov exponent Ai of system (|TT|) 
as a function of the damping parameter 7 and the frequency of modulation 
/i, for e = 0.01, F = 5 and Qp = u = 1 (the resonant case) and with the 
initial conditions: Rea = 10, Ima = 0. 
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Figure 9: Time dependence of Re a of system (ITT!) , for 7 = 0.01, /x = 0.15, 
F = 5, e = 0.01 and Qp = u = 1. The system starts from the initial 
condition: Rea = 10, Im a = 0. Chaotic beats. 



(nonhnear couplers) is described by the following hamiltonian: 

H = Ho + Hi + H2, (12) 

where: 



Ho 
Hi 
H2 



(13) 
(14) 
(15) 



The hamiltonian Hq represents the two single Kerr oscillators. Hi is the 
hamiltonian of the interaction between them and H2 describes the interac- 
tion of the two oscillators with the external fields. The values ai and 02 are 
complex dynamical variables, Uj denote the frequencies of the free vibrations 
of the two single oscillators; ei and €2 are Kerr parameters describing nonlin- 
earity in these sub-systems (oscillators) and 612 is the parameter of nonlinear 
coupling between them. 

The equations of motion for variables ai and 02 have the form: 



^^1 • • * 2 • * 

— — = — zwifli — leia^a-i — iei2a2aia2 + 

(JjL 



7i«i 



'•2 ■ • * 2 • * 

= —tU2a2 — «e2Ci2'3'2 ~ iei2aiaia2 + 



(16) 
(17) 



72^2 



where the terms 7iai and 720-25 describing the dissipation of energy, with 
damping constants 71 and 72, are added to the equations of motion on phe- 
nomenological grounds. In the autonomous and conservative case, that is 
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when 7i = 72 = Fi = F2 = the solutions of equations ( [T6|) - (fT7|) have the 
form: 



ai{t) 
a2(t) 



p~i{uJi+eia*f^aio+ei2a2Qa2o)t 



e 



-i(ui2+e2a2f)a20+i^l2a'loaio)t 



(18) 
(19) 



where ai(t) 



t=o- 



Oio and 02 (t) 



|(=o= 0.20 denote initial conditions. 
In the nonautonomous case (|T6l) -( fT71) we find periodic solutions in the form 
of the following functions (in analogy to the case of the single Kerr oscillator): 

-i{uii+eixlxi+ei2X2X2)t 



ai(t) 
a2{t) 



xie 
X2e 



-i{u}2+e2X2X2+ei2xlxi)t 



(20) 
(21) 



that is in the form of the solutions of the autonomous one. Functions fl20l) - 
( l2Ti) satisfy the equations of motion ( |T6l) -( |T7l) on condition that Xj = Fj/'jj, 
j = 1,2. As a result these solutions have the form: 



ai{t) 
02 (t) 



7i 

72 



exp 



exp 



/ f2 p2^ 

-i [ui + ei^ + ei2^ 
\ 7i 72 , 

-t\uj2 + 62— + ei2 — 
\ 72 I1 , 



t 



(22) 
(23) 



Then nip = 0Ji + eiFfhl + ei2Fihl and ^2^ = ^2 + e2Fihl + £12^2/72. 
In the phase plane (i?eai, Imai, Rea2, Ima2) the periodic solutions ( 122|) - 
satisfy the phase equations (limit cycles - circles): 



{Re ajY + {Im ajY - 
for any values of frequencies ui, 002- 



71' 



where j 



1,2, 



(24) 



3.2 The dynamics of the system in the phase space 

Coupling the single Kerr oscillator described by equation (jl]) to the other 
analogous oscillator but with different own frequency causes distinct changes 
in its dynamics, as illustrated by the diagrams in the phase space. 

Let us consider two nonlinearly coupled Kerr oscillators (sub-systems) 
described by equations f|T6l) - f|T7j) if ui = 1, uj2 = 0.5, ei = £2 = 0.01, 
Fi = F2 = 5, 71 = 72 = 0.5, fiip = 3.0, n2p = 2.5 and 612 = 0.01. Then, 
in compliance with fl2^ - fl2^ . the periodic solutions of equations ffTBl) - ffT7l) 
have the form: 

ai(t) = 10e-=^-°**, (25) 
a2{t) = lOe-2-5^*, (26) 

and in the phase space they satisfy the following equations: 

{Reaif + {Imaif = 100, (27) 
{Rea2y + ilma2f = 100. (28) 
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As a result, for the initial conditions aio = 10 and 020 = 10 the phase points of 
both subsystems draw the same circle of the radius r = 10 described by equa- 
tions (1271) - (!28l) . with frequencies Qip = 3.0 and = 2.5, respectively. But if 
the phase point representing the first subsystem (at frequency Ui = 1) starts 
from another initial condition instead of Oio = 10, and 020 is fixed (020 = 10), 
the following behaviour is observed: the phase point representing the first 
sub-system after some time tends to one of the two orbits being attractors of 
subsystem ([16]): {Reaif + {Imaif = 100 or (iJeai)^ + {Imaif = 6.6987. 
The second orbit belongs to the periodic solution: ai(t) = (10/(1 — i{2 — 
V3))exp{-3.0it) and a2(t) = {10/ {1 - i{2 - V^))exp{-2. bit). This situation 
is illustrated in Fig. UHl In Fig. ITOt^a) the phase point representing the 
first sub-system starts from the initial condition ReaiQ = 0, ImaiQ = 15 
(the initial condition for the second subsystem (1T7|1 is fixed: Re 020 = 10, 
Im a2o = 0) and after the time t = 50 it goes into the attractor of the radius 
r = 10, described by equation ( l27j) . However, FigJTOlfb) shows the phase 
point starting from the initial condition Re aio = 0, Im aiQ = 7 and going 
into the orbit of the radius r' = a/6.6987. 

The third solution of (16) and (17) is completely unstable and has the 
form: 

ai{t) = (10/(1 - i{2 + y/3))exp{-3.0it) 
and a2(t) = (10/(1 - i{2 + y^))exp{-2.5it). 

3.3 Basins of attraction 

FigHH shows two attractors and their basins of attraction of sub-system 
( lT6l) marked by appropriate colours; the initial condition of sub-system (1T71) 
is fixed (020 = 10). The yellow area marks the basin of attraction of the 
attractor of the radius r = 10, whereas the blue one corresponds to the 
basin of attraction of the attractor of the radius r' = a/6.6987. The basin 
marked in yellow has a special geometry: it consists of two separate areas, 
one of which has a spiral-like form. Both, the slip of the spiral and its width 
decrease towards moving away from the centre, similarly, as for the single 
Kerr oscillator. Contrary to the case of the single Kerr oscillator there is 
an island in the central part of the basin. The remaining area (blue colour) 
is the basin of attraction of the other attractor (the circle of the radius 
r' = V6.6987). 

The attractor of the radius r' = V6.6987 is stable (it is fully in its own 
basin of attraction), however the attractor of the radius r = 10 is semistable 
(it is partly in its own basin of attraction and partly in the basin of attraction 
of the other attractor). As a result the transition from the attractor of the 
radius r = 10 to the attractor of the radius r' = V6.6987 is possible, but that 
in opposite direction is impossible, because the phase point starting from any 
position on the circle {Re a)^ -|- {ImaY = 6.6987 always returns to it. 

The types of attractors change after changing the initial condition of 
the sub-system (1T71) . For example, FigjlH shows the attractors and their 
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Figure 10: The phase trajectories of sub-system (11611 if wi = 1, U2 = 0.5, 
61=62 = 0.01, Fi = Fa = 5, 71 = 72 = 0.5, fiip = 3.0, figp = 2.5, eia = 0.01 
and for the initial conditions: (a) Reaio = 0, ImaiQ = 15, -Re 020 = 10, 
Ima2o = 0; (b) Reaio = 0, Imaio = 7, Rea2o = 10, Ima2o = 0. 
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Figure 11: Basins of attraction of two coupled Kerr oscillators. The param- 
eters of the systems are: ui = 1, uj2 = 0.5, ei = €2 = 0.01, Fi = F2 = 5, 
7i = 72 = 0.5, flip = 3.0, Q2p = 2.5, = 0.01, and the initial conditions 
are: ReaiQ = 0, ImaiQ = 15, Rea2o = 10, Ima2o = 0. Stable (r' = 1/6. 6987) 
and semistable (r = 10) attractors. 



basins of attractions of sub-system f|T6l) for a2o = 0. The basin of attraction 
corresponding to the circle of the radius r = 10 is marked by yellow; the 
remaining area (blue colour) refers to the basin of attraction of the circle 
of r' = -\/6.6987. As we can see, in this case the attractor with radii (r' = 
Vs. 6987) is stable, and the other one is completely unstable. 

3.4 Generation of chaotic beats 

Globally, the behaviour of system (IT6l)-(IT71) is presented in FigJTSTa) show- 
ing the Lyapunov map in the parameters space (€12,7). We find that strong 
chaotic behaviour of the system is much common that for the single Kerr 
system. We also notice that if we increase the dumping in the system we 
must increase the coupling between subsystems to achieve the chaotic be- 
haviour. The full spectrum of the Lyapunov exponents {Ai, A2, A3, A4} versus 
612 shows the regions of order or chaos in the cross section of the map for 
the damping parameter 7 = 71 = 72 = 0.5 (FigJT3](b)). If Ai > then the 
system is chaotic, and if Ai < 0, it behaves periodically. The system with the 
parameters of Fig. [13] and for the coupling constant €±2 > 1.6 manifests ex- 
tremely unstable behaviour and its solutions are divergent to infinity. There 
is also a region of hyperchaotic behaviour of the system in which two highest 
Lyapunov exponents are positive (for ei2 > 0.44). 

Taking from Fig{T3t^a)-(b) the appropriate values of the damping constant 
7i = 72 = 0.5 and the nonlinear coupling between the Kerr oscillators ei2 = 
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Figure 12: Basins of attraction of two coupled Kerr oscillators. The pa- 
rameters are the same as in Fig. [Tl|, but Rea2o = 0, Ima2Q = 0. Stable 
(r' = -\/6.6987) and unstable (r = 10) attractors. 



0. 45 we can generate chaotic beats in the system of two nonlinearly coupled 
Kerr oscillators f lT6|) - f|T7|l . both being initially in the periodic state. Such 
chaotic beats are shown in FigJH] illustrating the time dependence of Re ai 
and Re 02- Because the spectrum of Lyapunov exponents of the beats shown 
in FigHH is 

{0.0708, 0.0205, -1.4552, -1.5215} 

and contains two positive values we can even call it hyperchaotic beats 

m- 

4 Conclusion 

The Kerr effect and the Kerr couplers considered in this paper have great 
potential in studies and applications of optical devices like optical fibres 
or couplers. From this point of view the main results of this paper are: 

1. Tunnelling properties of periodic attractors (mutual inter-penetration of 
attractors and basins of attraction) lead to the possibility of switching of Kerr 
oscillator system between different semi-stable attractors by changing initial 
condition (Sec. 2. 3). 2. As shown in Sec. 2. 4 the dynamics of the Kerr system 
can be controlled by appropriate switching of the parameters of the system or 
the external field. Temporary changes in parameters also switch the system 
between different periodic states. 3. The system with Kerr nonlinearity 
is able to generate chaotic beats (Sec. 2. 5). The properties of these beats 
depend on specific choices of parameters of the system and the external field. 
4. For two coupled Kerr oscillators it is possible to change the stability of 
attractors of a given sub-system by changing the initial conditions of the 
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Figure 13: (a) The Lyapunov map of the system ( fT6l) -( fT7|) in the parametric 
space (ei2,7) where 7 = 71 = 72- (b) The full spectrum of Lyapunov expo- 
nents {Ai, A2, A3, A4} of the system (|T6l)-(|T71) as a function of the nonlinear 
coupling parameter eu and for 7 = 71 = 72 = 0.5. The rest of parameters 
for both figures are: Ui = 1, W2 = 0.5, ei = €2 = 0.01, Fi = F2 = 5, Qip = 3 
and fl2p = 2.5. The system starts from the initial conditions: Re ai = 10, 
Im ai = 0, Re 02 = 10, Ima2 = 
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Figure 14: Time dependence of Re ai (a) and Re a2 (b) of the system (IT6l)- 
f lT7|) . Parameters are the same as in Fig JTST b). and ei2 = 0.45. Chaotic 
beats. 
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other sub-system and the dynamics of one sub-system can be controlled by 
changing the initial conditions of the other one (Sec.3.3). 5. Moreover, 
chaotic (or hyperchaotic) beats in the system of two coupled Kerr oscillators 
can be generated (Sec. 3.4). 
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